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ABSTRACT 

o . 

' Context. The mass-loss mechanism in normal K-M giant stars with small variability amplitudes is not yet understood, although they 

are the majority among red giant stars. 

in : 

Aims. We present high-spatial and high-spectral resolution observations of the M7 giant BK Vir with AMBER at the Very Large 
Telescope Interferometer (VLTI). Our aim is to probe the physical properties of the outer atmosphere by spatially resolving the star 
in the individual CO first overtone lines. 

Methods. BK Vir was observed between 2.26 and 2.31 /jm using the 16-32-48 m telescope configuration with an angular resolution 

Mh of 9.8 mas and a spectral resolution of 12000. 

Results. The uniform-disk diameters observed in the CO first overtone lines are 12-31% larger than those measured in the continuum. 
(— I ' We also detected asymmetry in the CO line-forming region, which manifests itself as non-zero/non-7r differential and closure phases. 

The data taken 1.5 months apart show possible time variation on a spatial scale of 30 mas (corresponding to 3 x stellar diameter) at 
the CO band head. Comparison of the observed data with the MARCS photospheric model shows that whereas the observed CO line 
spectrum can be well reproduced by the model, the angular sizes observed in the CO lines are much larger than predicted by the model. 
Our model with two additional CO layers above the MARCS photosphere reproduces the observed spectrum and interferometric data 
in the CO lines simultaneously. This model suggests that the inner CO layer at ~1.2 R+ is very dense and warm with a CO column 
density of ~10 22 cm" 2 and temperatures of 1900-2100 K, while the outer CO layer at 2.5-3.0 R* is characterized by column densities 

'— 1 ■ of 10 19 -10 20 cm" 2 and temperatures of 1500-2100 K. 

Conclusions. Our AMBER observations of BK Vir have spatially resolved the extended molecular outer atmosphere of a normal M 
giant in the individual CO lines for the first time. The temperatures derived for the CO layers are higher than or equal to the uppermost 
layer of the MARCS photospheric model, implying the operation of some heating mechanism in the outer atmosphere. This study 
also illustrates that testing photospheric models only with the spectra of strong molecular or atomic features can be misleading. 
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t— ( | 1. Introduction 1987 stars classified as semi-regular or irregular giants (denoted 

y—i _ , , . „ , , as SRA, SRB, SRD, LB, and L in the catalog), in contrast to 

j. . The mass loss in late evolutionary stages of low- and 1524 starf , classified as Mira variables (see , however , Lebz elter 

e > ■ intermediate-mass stai-s (1 £ M/Mq <= 8), particularly in the red gt a , 1995 for ^ uncertainty in the c i assiflcat i on m the GCVS). 

_ giant branch (RGB) and asymptotic giant branch (AGB) is im- There ^ ^ sem i-regular or irregular variables whose 

■ portant n0t 0nly fOT the evolutlon of the star ltself but also for the spectral type is unknown in the catalog. Therefore, normal M gi- 
rt • chemical enrichment of the interstellar medium. The mass loss antfj wkh smaU variabilit amplitudes may outnumb er Mira stars 

in Mira-type stars in the AGB phase, which are characterized considerably 

by the large-amplitude pulsation with a variability amplitude of n . • „ • ,.. ■ ■ , f 

, _ & F , ■ , r , , t ■ Despite their small variability amplitudes and the absence or 

AV * 6 mag and a period of ~1 year, has been observation- ^ erio dicities, normal M iants are experiencing mass loss 

ally and theoretically intensively studied. However, normal (i.e., .., „ , . 7, . • , .■ „ u . 

3 ... . J . , J ..... .. . . .T, with mass-loss rates comparable to those in (optically bright) 

non-Mira-type) M giants show variability amplitudes or AV = , ,■ . v T , „ ■ , • ., A/ ^ D , , 

, _ { „ , • , , Mira stars. Normal M giants in the AGB phase show mass- 

1-2 mag, much smaller than those of Mira stars, without clear ^ m ^ of 10 -7_ 10 -6 M -i and expansion velocities of 5- 

penodic variations. They are classified as semi-regular ^or irreg- 15 km ^ ( Gonzalez Delgado et al. 2003; Winters et al. 

ular variables. Among the M stars listed in the General Catalog i onn o i t-, r> i * i hmni -m. • 1 e ■ -c 

„ A ~ ,a , r^T^i ■ 2003, De Beck et al. 2010 . There is also evidence or significant 

of Variable Stars (GCVS ) 4.2 (Samus et al. EES, there are ^, oss ^ RGR starg ^ mass . loss rates up tQ 1Q _ 6 ^ yr _! 

Send offprint requests to: K. Ohnaka near the li P of the RGB (° ri g lia et al-MMl Ita et al.|2SD. 

* Based on AMBER observations made with the Very Large The mass-loss rate (M) and expansion velocity (t, exp ) of our 

Telescope Interferometer of the European Southern Observatory, target, BK Vir, are estimated to be (1.5 - 4) X 10 7 M q yr 1 

Program ID: 081.D-0233(A) (AMBER Guaranteed Time Observation) and 4-7.5 km s~', respectively (Gonzalez Delgado et al. 2003, 
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Table 1. Summary of the AMBER observations of BK Vir with the E0-G0-H0 AT configuration. 

# Night f obs Bp PA Seeing Airmass DIT Number Calibrator 
(UTC) (m) (°) Q (s) of frames 

1 2009 Jan 19 07:43:32 13.9/27.9/41.9 73 0.9 1.2 6 75 vHya T 

2 2009 Feb 28 06:15:16 15.6/31.2/46.8 74 0.7 1.1 6 75 vHya+ 

3 2009 Apr 15 03:43:30 15.9/31.8/47.8 74 0.8 1.1 6 75 BCrv 

Note: r means that the calibrator was observed before and after the science target. Seeing is in the visible. 



Winters et al. 2003 ), which are approximately the same as in the 
prototypical Mira o Cet with M = (1 - 8) X 10~ 7 Mq yr^ 1 and 
^exp = 1-5-7 km s _1 (Winters et al. 2003 ). This illustrates the im- 
portance of the contribution of normal M giants to the chemical 
enrichment of the interstellar medium. 

For understanding the mass-loss mechanism in red giant 
stars, it is indispensable to probe the physical properties of the 
region between the upper photosphere and the innermost part 
of the circumstellar envelope, where the energy and momen- 
tum are expected to be deposited for the wind acceleration. It 
is now known from spectroscopy and high-spatial resolution ob- 
servations in the IR that there is a dense, warm (~ 1000-2000 K) 
molecular outer atmosphere, the so-called MOLsphere, extend- 
ing to a few R+ not only in Mira stars but also in normal K-M 
giants with small variability amplitudes (e.g., Tsuji et al. 119971 
Tsuji 200T1 P errin e t al. l20U4l Ohnaka|20W! Ohnaka et al. l20031 
Takami et al. 2009). This MOLsphere is considered to play an 
important role in driving mass outflows. In the case of Mira stars, 
the large-amplitude pulsation can fairly explain the presence of 
the extended atmosphere (e.g., Ohnaka et al. 2006; Wittkowski et 
al. 120071 12008; Woodruff et al. 120091 However, the amplitudes 
of the stellar pulsation are much smaller in normal K-M giants, 
and therefore, the origin of the MOLsphere in these stars is by 
no means clear. In addition to the MOLsphere, the detection of 
UV emission lines as well as the Ha line suggests the presence 
of a chromosphere (e.g., Eaton [1995l Isabel Perez Martinez et al. 
1201 H and references therein). This means that there are both hot 
and warm components in the outer atmosphere of normal K-M 
giants, and the spectral analyses of the CO lines in the UV and 
IR support this picture (e.g., Wiedemann et al. 1994; McMurry 
& Jordan[2000]). 

The multi-component nature of the outer atmosphere is 
a key to understanding the mass-loss mechanism in red gi- 
ants. IR spectro-interferometry, which enables spatially resolved 
spectroscopy of molecular or atomic spectral features with a 
spatial resolution much higher than achieved with single-dish 
telescopes, is effective for probing the physical properties of 
the outer atmosphere. While Mira stars have been intensively 
studied by IR spectro-interferometry as mentioned above, IR 
spectro-interferometric observations of normal red giants are 
still very scarce. Recently, Marti- Vidal et al. (1201 II ) have pre- 
sented near-IR spectro-interferometric observations of the M6/7 
giant RS Cap in the CO first overtone bands and H2O bands with 
a medium spectral resolution of 1500 using VLTI/ AMBER. In 
this paper, we report on VLTI/ AMBER observations of the M7 
giant BK Vir in the individual CO first overtone lines near 2.3 fim 
with a higher spectral resolution of 12000. BK Vir shows no pe- 
culiar features and, therefore, is a good representative of normal 
M giants (its evolutionary status is discussed in Sect. 14. U . It was 
selected also for technical reasons, such as the brightness and 
angular size appropriate for AMBER observations. 

The paper is organized as follows. Our AMBER observa- 
tions and data reduction are described in Sect. [2] The observa- 
tional results are shown in Sect. [3] followed by the modeling of 
the data presented in Sect. [4] The implications of the results of 



the modeling are discussed in Sect. [5] Conclusions are presented 
in Sect. |6] 

2. Observations 

AMBER (Petrov et al. [20071 ) is the near-IR (1.3—2.4 /j.m) 
spectro-interferometric instrument at VLTI. It enables us to ob- 
tain information on the amplitude and phase of the complex 
Fourier transform of the object's intensity distribution with a 
spectral resolution of up to 12000 by combining three 8.2 m Unit 
Telescopes (UTs) or 1.8 m Auxiliary Telescopes (ATs). AMBER 
measures several quantities related to the complex Fourier trans- 
form of the object's intensity distribution: visibility, differen- 
tial phase (DP), and closure phase (CP). The visibility (or also 
called visibility amplitude) contains information about the size 
and shape of the object. The DP can measure the deviation of the 
object's phase in a spectral feature from that in the continuum, 
and non-zero DP within a spectral feature means that the spectral 
feature-forming region is asymmetric. The CP is the sum of the 
measured phases around a closed triangle of baselines. The CP 
is always zero or n for point-symmetric objects, and non-zero 
and non-7r CPs indicate an asymmetry in the object. 

BK Vir was observed at 3 epochs between January and April 
2009 with AMBER using the E0-G0-H0 telescope configuration 
with ATs, as summarized in Table [TJ We used the high-spectral 
resolution mode in the K band (HR_K) with a spectral resolution 
of 12000 in the wavelength range between 2.26 and 2.31 //m, 
which includes the CO (2,0) band head at 2.294 fim. The VLTI 
fringe-tracker FINITO was used to stabilize the fringes for a 
Detector Integration Time (DIT) of 6 s. This long DIT was nec- 
essary to achieve reasonable SNRs in the CO lines, where the 
correlated flux (i.e., flux x visibility amplitude) is very low. 

We reduced the data with amdlib ver2._fl. We took the best 
80% of the frames in terms of the fringe SNR (Tatulli et al. 
120071 . The interferometric calibrators (v Hya and /3 Crv) were 
observed just before and/or after BK Vir. We adopted angular 
diameters of 3.40 and 4.45 mas for v Hya and (3 Crv, respec- 
tively (Richichi & Perschron 2005). The errors of the calibrated 
visibilities, DPs, and CPs were estimated as described in Ohnaka 
et al. (120091 . 

The AMBER data taken in HR_K mode before December 
2009 are affected by the Fabry-Perot effect originating in the 
dichroics of the InfraRed Image Sensor (IRIS), which stabilizes 
the image motionQ. This effect appears as time-dependent, high- 
frequency beating in raw visibilities, DPs, and CPs, when plotted 
as a function of wavelength. Fortunately, this effect is mostly re- 
moved by using the calibrators observed just before and/or after 
the science target, although the residual of this effect appears as 
the small-amplitude beating in some cases (as seen in the ob- 
served CP at wavelengths below 2.293 /mi in Fig. QJ). 

The long DIT of 6 s makes the absolute calibration of the ob- 
served visibilities uncertain, because we cannot know the frac- 

1 Available at http://www.jmmc.fr/data_processing_amber.htm 

2 AMBER data taken after the replacement of the IRIS dichroics in 
December 2009 do not show this effect any more. 
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tion of time during DIT when the fringes were lost, and the 
fringe tracking may work differently for the science target and 
the calibrator. Therefore, we took the following approach for the 
absolute calibration of the observed visibilities. We note that the 
TT-band uniform-disk diameters of BK Vir reported in the liter- 
ature show little time variation within the measurement errors: 
11.2 + 0.6 mas (Dyck et al. [19981 . 10.73 + 0.2 3 mas (Perrin 
et al. [19981 , and 11.49 + 0.40 mas (Perrin et al. 120031 . In the 
present work, we adopted a uniform-disk diameter of 10.73 mas, 
which has the smallest error among the above measurements. We 
scaled the observed calibrated visibilities in the continuum to 
those expected from the uniform disk with an angular diameter 
of 10.73 mas. 

The wavelength calibration was carried out using the telluric 
lines identified in the spectra of the calibrators. As a template 
of the telluric lines, we convolved the atmospheric transmis- 
sion spectrum measured at the Kitt Peak National Observatory^ 
to match the spectral resolution of AMBER. The uncertainty 
in wavelength calibration is 4 x 10~ 5 //m (~5.2 kms _I ). The 
wavelength scale was converted to the laboratory frame using 
a heliocentric velocity of 16.5 km s _1 measured for BK Vir 
(Gontcharov l2006l Famaey et al. 120051 . 

Since no spectroscopic standard star was observed, we used 
the interferometric calibrators for the spectroscopic calibration 
to remove the telluric lines from the observed spectrum of 
BK Vir. However, both /3 Crv (G5II) and v Hya (K0/K1III) show 
CO absorption lines. In this case, the calibrated spectrum of 
the science target is derived as Ff 1 = F^g /(F™£ / f^ al ), where 
F sci(cai) d F sci(cai) denote the tm and b se rved (i.e., includ- 

* obs v ' 

ing the atmospheric transmission and the detector's response) 
spectra of the science target (or calibrator), respectively. To esti- 
mate the true spectrum of v Hya, we used the high-resolution K- 
band spectrum of a Boo (K1III) obtained by Wallace & Hinkle 
(1996), which has a similar spectral type to v Hya. The spec- 
trum of a Boo, which has a spectral resolution of 89000, was 
convolved down to match the resolution of AMBER. This con- 
volved spectrum of a Boo reproduces the strength of the CO 
band head observed in v Hya, indicating that it can be used as 
the true spectrum of v Hya. 

For fi Crv, however, no high-resolution CO spectra of 
stars with similar spectral types are available in the litera- 
ture. Therefore, we estimated the true spectrum for ft Crv by 
constructing a model CO spectrum using the MARCS models 
(Gustafsson et al. 2008). The MARCS models represent spher- 
ically symmetric photospheres with molecular and atomic line 
opacities taken into account using the opacity sampling method. 
Each MARCS model is specified by effective temperature (T e ff), 
surface gravity (g), microturbulent velocity (UmkroX chemical 
composition, and stellar mass (M*). These parameters derived 
for /? Crv are T e s - 5 100-5 145 K, log g = 2.2-2.7 (given in units 
of cm s~ 2 throughout the paper), i> m i cl0 — 1.5-2.1 kms^ 1 , and 
3.0- 3.7 Mq (Luck et al. [1^951 Takeda et al. 120081 Lyub imkov 
et al. l20T0l i. The CNO abundances derived by Luck et al. ([19951 
suggest the mixing of CN-cycled material. We selected a model 
atmosphere with parameters as close to these values as possible 
from the MARCS model gricQ. The selected model is character- 
ized with r e ff = 5250 K, \ogg = 2.5, u m i C i- = 2 km s _1 , M+ = 
1 Mq, and "moderately CN-cycled" chemical composition with 
C/N = 1.5. The stellar mass of this model is smaller than the ob- 
servationally estimated values above. However, the atmosphere 

3 http://www.eso.org/sci/facilities/paranal/instruments/isaac/tools/ 
spectra/atmos_S _K.fi ts 

4 http://marcs.astro.uu.se 



is approximately plane-parallel with a geometrical thickness of 
~1% of the stellar radius. In this case, the stellar mass does not 
affect the atmospheric structure noticeably. Using the temper- 
ature and pressure distributions of this model, we computed a 
synthetic spectrum for the CO first overtone lines using the line 
list of Goorvitch (1994). The carbon abundance was adjusted 
to fit the CO band head in the uncalibrated spectrum of ft Crv. 
With F'jf constructed in this fashion, the calibrated spectrum 
of BK Vir was derived using the above equation. Although the 
strengths of the CO lines in two calibrators are very different (the 
line depths are 20-30% and 10% for v Hya and ft Crv, respec- 
tively), the spectra of BK Vir calibrated with these stars agree 
well (see Fig. [3^, where both spectra calibrated with v Hya (2009 
Feb 28) and ft Crv (2009 Apr 15) are shown). This indicates that 
the calibrated CO line spectra of BK Vir are not significantly 
affected by the above, unusual spectroscopic calibration. 

3. Results 

Figure[T]shows the observed visibilities, DPs, CPs, and spectrum 
from the data set #3. The results from the other data sets are sim- 
ilar, and therefore, not shown. The figure reveals significant sig- 
natures of the CO lines in the observed visibilities, DPs, and CPs. 
Thanks to the long DIT of 6 s, the visibilities as low as 0.01 could 
be measured in the CO lines. As plotted in Fig.QJl, the uniform- 
disk diameter obtained from the observed visibilities shows that 
the star appears larger in the CO lines (12-14 mas) than in the 
continuum (10.73 mas) by 12-31%. However, the uniform-disk 
fit in the CO lines is poor with reduced x 2 values much higher 
than 1 (as high as 10 2 -10 3 ), which means that the uniform disk 
is a poor representation of the object's shape in the CO lines. 
Therefore, the uniform-disk diameter in the CO lines should be 
regarded as a very rough estimate of the characteristic size of the 
object. 

Figures QJ-[TJi show the detection of non-zero DPs (partic- 
ularly clear on the 47.8 m baseline) and non-zero/non-7r CPs 
in the CO lines, which indicate an asymmetry in the CO line- 
forming region. High-spatial resolution observations have re- 
vealed that asymmetries in the atmosphere of red giant stars 
are common. The images of the Mira star R Cas in the TiO 
bands obtained by Weigelt et al. (1996) show elongation, while 
the image in the continuum appears circular. The optical im- 
ages of the prototypical Mira star o Cet taken with the Hubble 
Space Telescope show significant asymmetries (Karovska et al. 
1997). Tuthill et al. (119991 ) also revealed departures from circu- 
lar symmetry in five Mira stars (including the above two stars) 
in the optical and near-IR. The aperture synthesis images of 
the Mira stars R Aqr and U Ori in the H band obtained with 
the Infrared Optical Telescope Array (IOTA) show asymmet- 
ric, extended atmospheres (Ragland et al. 2008; Pluzhnik et al. 
2009 ). U Ori shows asymmetry in the K band as well (Mondal & 
Chandrasekhar 2004). There are more stars for which asymme- 
tries are indicated from the detection of non-zero CPs. Bedding 
et al. (1997 ) detected non-zero CPs for the Mira-like star R Dor. 
Ragland et al. (2006) also report asymmetries in a sample of 
AGB stars based on the detection of non-zero CPs. More re- 
cently, Wittkowski et al. (2011) have shown asymmetries in the 
CO first overtone bands in Mira stars by AMBER observations 
with a spectral resolution of 1500. Our AMBER observations 
are the first study to detect asymmetry in a red giant star in the 
individual CO lines. 

Figure^, which shows an enlarged view of the visibilities in 
four representative CO lines obtained on the 47.8 m baseline, re- 
veals that the observed visibilities are approximately symmetric 
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Fig. 1. AMBER data of BK Vir (data set #3). In each panel, the scaled observed spectrum is plotted by the black solid lines. a)-c) 
Visibilities observed on the 15.9, 31.8, and 47.8 m baselines (red lines). The corresponding spatial resolutions are also given, d) 
Uniform-disk diameter (red line) derived by fitting the visibilities shown in the panels a)-c). e)-g) Differential phases observed on 
the 15.9, 31.8, and 47.8 m baselines (red lines), h) Closure phase (red line). 



with respect to the line center. On the other hand, Fig.fJJ) shows 
the visibility of the red supergiant Betelgeuse observed on the 
1 1 .54 m baseline (data set #5 of Ohnaka et al. 1201 11 1. Because the 
angular diameter of Betelgeuse (42.5 mas, Ohnaka et al. 1201 11 1 
is four times larger than that of BK Vir (10.73 mas), the 1 1.54 m 
baseline for Betelgeuse probes nearly the same spatial scale in 
terms of the stellar angular diameter as the 47.8 m baseline for 
BK Vir (0.9 x stellar angular diameter). Figure[2|3 clearly shows 
that the observed visibilities are asymmetric with respect to the 
line center in the case of Betelgeuse: the visibilities are charac- 
terized by the minima in the blue wing and the maxima in the red 
wing. These asymmetric visibilities are interpreted as temporally 
variable, inhomogeneous gas motions within 1.3-1.47?* over the 
stellar surface with velocity amplitudes of up to 20-30 km s _1 
(Ohnaka et al. 2009] 1201 11 1. The absence of such asymmetric 
visibilities in BK Vir suggests that the amplitude of the inhomo- 
geneous velocity field is much smaller than the AMBER's spec- 
tral resolution of 12000, say, ~5 km s _1 , and/or the spatial scale 
of the inhomogeneous gas motions is much smaller than the spa- 
tial resolution of the current AMBER data (9.8 mas). However, 
the former possibility is more favorable, because the non-zero 
DPs and non-zero/non-7r CPs shown in Fig. Q] indicate that the 
inhomogeneity is spatially resolved with the current spatial res- 
olution. 



We examined time variation using the AMBER data taken at 
nearly the same uv points. The shortest and middle baselines of 
the data sets taken on 2009 Feb 28 and 2009 Apr 15 are almost 
the same (the 1 m difference in the longest baseline is too sig- 
nificant to study time variations). A comparison between these 
two data sets is shown in Fig. [3] The sampled wavelengths are 
slightly different for the two nights because of the difference in 
the correction to convert the observed wavelength scale to the he- 
liocentric scale. Figures [3^ and[3j; show that the observed spec- 
tra and the visibilities on the 32 m baseline measured on the 
two nights form continuous curves, suggesting no time varia- 
tion. This is also the case for the 16 m visibility shown in Fig.[3j5 
except for the CO band head at 2.2936 yum, where the visibil- 
ity measured on Apr 15 is 12-20% lower than that measured on 
Feb 28. This suggests possible time variation on a spatial scale 
of 30 mas (3 x stellar diameter). The decrease in the visibility 
corresponds to an increase in the uniform-disk diameter by 22% 
(from 16 to 19.5 mas) or an increase in the Gaussian FWHM 
by 14% (from 10.5 to 12 mas) between Feb 28 and Apr 15. 
However, the difference is within 3-cr of the data from Feb 28. 
Therefore, more observations are needed to definitively confirm 
time variations in the visibility. 
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Fig. 2. Visibilities in four representative CO first overtone lines 
in the red giant BK Vir (panel a) and the red supergiant 
Betelgeuse (panel b). In both panels, the black solid lines repre- 
sent the scaled observed spectra. The red solid lines represent the 
visibility obtained for BK Vir on the 47.8 m baseline (panel a) 
and that obtained for Betelgeuse on the 1 1.54 m baseline (panel 
b). 



4. Modeling of the AMBER data 

4.1. Determination of stellar parameters 

To interpret the AMBER observations, we use the MARCS mod- 
els. As mentioned in Sect. [2 these models assume spherical 
symmetry. Whereas the asymmetry in the CO-line-forming re- 
gion is detected, we use these spherical models as the first ap- 
proximation. To select a MARCS model appropriate for BK Vir, 
we determined its stellar parameters (T^, g, u m i cro , M+, and 
chemical composition) as follows. 

We estimated the effective temperature from the observed 
angular diameter and the bolometric flux obtained by integrat- 
ing the photometric data available in the literature (Two-Micron 
Sky Survey, Ne ugebau er & Leighton[l969 l NO MAD Catalog, 
Zacha rias et al. l2004l 2MASS, Sk rutskie et al. l2006l Fouque 
et al. IT9921 Kerschbaum & Hron [19941 IRAS Point Source 
Catalog). The photometric data were corrected for the interstel- 
lar extinction using E(B -V) = 0.023 (Schlegel et al. [T998l > and 
assuming Ay = 3.1E(B - V). Combining the angular diameter 
of 10.73 mas adopted in Sect.|2]and a derived bolometric flux of 
2.80 xl0~ 9 WirT 2 results in an effective temperature of 2920 K. 
While this agrees with the effective temperature of 3074 + 141 K 
derived by Dyck et al. (1998), they obtained a higher bolometric 
flux of 3.9 x 10~ 9 WirT 2 . If we adopt, as the uncertainty in the 
bolometric flux, a half of the difference between the values from 
Dyck et al. ( 1998) and from the present work, the total error (i.e., 
error resulting from the uncertainties in the angular diameter and 
in the bolometric flux) in our effective temperature is +150 K. 
Using the distance of 180 pc based on the Hipparcos parallax of 
5.53 + 0.68 mas (van Leeuwen 2007) and the above bolometric 
flux of 2.80 x 10~ 9 WirT 2 , we derived a luminosity of 2700 L© 
(Mbql = -3.8). With the error in the parallax and the above error 
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Fig. 3. Spectra (panel a) and visibilities (panels b: 16 m baseline 
and c: 32 m baseline) obtained on 2009 Feb 28 and 2009 Apr 15. 
In all panels, the data taken on 2009 Feb 28 and 2009 Apr 15 are 
plotted by the black and red dots, respectively. 



in the bolometric flux, the total uncertainty in the luminosity is 
+850 L e . 

To estimate the stellar mass, we compared the position of 
BK Vir on the H-R diagram with theoretical evolutionary tracks. 
Figure|4]shows evolutionary tracks for 1 and 2 Mq stars (Herwig 
20033 an d Bertelli et al. 120081 respectively), together with the 
observationally derived position of BK Vir. The figure suggests 
that the mass of BK Vir is close to 1 Mq. The adoption of a 
stellar mass of 1 Mq translates into a surface gravity of log g = 
-0.17 with a stellar radius of 201 Rq (0.94 AU). 

Whereas BK Vir is classified as an AGB star in most of the 
literature, Fig. [4] shows that the luminosity of BK Vir is just at 
the tip of the RGB. Moreover, Lebzelter & Hron (1 19991 1 report 
negative detection of 99 Tc, which would be strong evidence of 
the third dredge-up in the AGB. Therefore, we cannot conclude 
whether BK Vir is an RGB star or an early- AGB star that has not 
yet experienced the third dredge-up. In either case, however, the 
surface chemical composition is expected to be changed by the 
first dredge-up in the RGB phase, which mixes the CN-cycled 
material to the stellar surface. This makes the MARCS models 
with the "moderately CN-cycled" chemical composition (with 
[Fe/H] = 0.0) suitable for BK Vir. 

No analysis of the micro-turbulent velocity in BK Vir is 
available in the literature. However, the high-resolution spec- 
troscopic analysis of M7-8 giants with T e ff 5 3000 K by 
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Fig. 4. Observationally derived position of BK Vir on the H-R 
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AGB stage (early AGB and thermally-pulsing AGB) is marked 
in blue and red for the evolutionary tracks for 1 Mq and 2 Mq, 
respectively. 



We 



Tsuji ( 2008 ) shows micro-turbulent velocities of 3-4 km s 1 . 
adopt these values for BK Vir as well. 

4.2. Comparison with the MARCS photospheric model 

The MARCS model with the parameters as close to the derived 
values as possible is characterized by T e g = 3000 K, logg = 
0.0, M+ = 1 Mq, Umic-o = 2 km s _1 (models with i> m i cro = 3- 
4 km s are not available for this combination of T e g, log g, M+, 
and chemical composition, and the model with u m i cro = 2 km s _I 
is the closest one). Using the temperature and pressure distri- 
butions of the downloaded MARCS model, we computed the 
monochromatic intensity profile and spectrum over the observed 
wavelength range using the CO line list of Goorvitch ( 1994 ), as 
described in Ohnaka et al. (120061 ) . The monochromatic visibility 
was calculated by the Fourier transform of the monochromatic 
intensity profile. The intensity profile, visibility, and spectrum 
were then convolved with the AMBER's spectral resolution. The 
angular scale of the model was determined so that the uniform- 
disk diameter derived from the model visibilities in the contin- 
uum matches the adopted value of 10.73 mas. 

Figure [5] shows a comparison between the observed data 
(data set #3) and the predictions by the MARCS photospheric 
model. While the observed CO spectrum is well reproduced, the 
predicted visibilities in the CO lines are too high compared to 
the observed data. We also used other available MARCS models 
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Fig. 5. Comparison between the observed interferometric data 
(data set #3) and the MARCS model. In all panels, the black 
thin lines represent the predictions of the MARCS model with 
the parameters most appropriate for BK Vir (see Sect. 14.21 ). The 
thick colored lines represent the observed data, a) Spectrum. 
Residuals of the strong telluric lines are marked by ffi. b)-d) 
Visibilities on the 15.9, 31.8, and 47.8 m baselines, respectively, 
e) Closure phase. 



with slightly different stellar parameters (T e ff = 2700-3100 K, 
log g = -0.5 - +0.5, M* = 1-2 Mq, and u micro = 2-5 km s -1 ), 
but they cannot reproduce the observed visibilities in the CO 
lines, either. Figure [6] shows the intensity profiles and visibili- 
ties at three representative wavelengths in the continuum, at the 
CO band head, and in an isolated line. The figure illustrates that 
the photosphere extends to 1.17 R* in the CO band head and in 
the isolated line, but this geometrical extension is insufficient to 
explain the observed visibilities (Fig. [6J;). The too high visibili- 
ties in the CO lines predicted by the MARCS model means the 
actual CO-line-forming layer is much more extended than the 
model predicts or there is an additional component contribut- 
ing to the CO lines above the photosphere modeled by MARCS. 
Therefore, our AMBER observations have spatially resolved the 
MOLsphere in a normal M giant for the first time in the individ- 
ual CO lines. The same trend was reported for the normal M6-7 
giant RS Cap by Marti- Vidal et al. d2^TT1) based on AMBER ob- 
servations of the CO bands with a lower spectral resolution of 
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1500. These studies illustrate that testing photospheric models 
only with the spectra of strong molecular or atomic features can 
be misleading. 

4.3. Modeling of the MOLsphere 

To characterize the physical properties of the MOLsphere, we 
interpreted the AMBER data with a model in which one or two 
extended CO layers with constant temperatures and densities are 
added above the MARCS photospheric model. This model is the 
same as that described in Ohnaka (2004), but with the follow- 
ing modification: each layer is assumed to be geometrically thin 
with a thickness of 0.1 R+. The inner radius of each layer is 
used as the parameter to represent the size of the layer. The ra- 
dius, temperature, and column density of each layer are changed 
as free parameters to fit the observed spectrum and visibilities 
simultaneously. We tentatively adopted a microturbulent veloc- 
ity of 4 kms~' for the CO layers based on the spectroscopic 
analysis of Tsuji (2008) for M7-8 giants. We first attempted to 
explain the observed data with models with one additional CO 
layer but could not reproduce the spectrum and visibilities si- 



multaneously. Therefore, we attempted to explain the data with 
the models with two CO layers. Table [2] gives the range of the 
parameters searched in our modeling. Note that the radius of the 
inner layer was set to be equal to or larger than 1.2 /?* to avoid 
an overlap with the MARCS photospheric model, which extends 
to 1.17/?*. 

Figure [7] shows a comparison between the observed data 
(data set #3) and the best-fit model. The model intensity pro- 
files and visibilities at three representative wavelengths (contin- 
uum, CO band head, and isolated CO line) are plotted in Fig. [8] 
This model is characterized by the inner CO layer at 1.2 /?* 
with 2000 K and a CO column density of 2 x 10 22 crrr 2 and 
the outer CO layer at 2.5 /?* with 1900 K and a column density 
of 2 x 10 19 cm 2 . The column densities of the inner and outer CO 
layer correspond to 16% and 0.016% of the CO column density 
in the MARCS photospheric model. We found the range of the 
radius, temperature and CO column density of the inner layer 
to be 1.2-1.25 /?*, 1900-2100 K, and (l-2)xl0 22 cirT 2 , respec- 
tively. The derived ranges for the radius, temperature, and CO 
column density of the outer layer are 2.5-3.0 /?+, 1500-2100 K, 
and 10 19 -10 2() cm -2 , respectively. As shown in Fig. [7] the ob- 
served visibilities and spectrum are well reproduced. The model 
also predicts the CP to jump in the CO lines, in particular at 
the CO band head, as seen in the AMBER data. This can be 
explained as follows. Figure [SJ) shows that the inner CO layer, 
which is optically thick, makes the star appear larger in the CO 
band head than in the continuum by 30%. Then the longest base- 
line samples the second visibility lobe in the CO band head, 
as shown in Fig. [8J; (red line and red triangle at a spatial fre- 
quency of ~100 arcsec -1 ). The phase in the second visibility 
lobe is 7r, while the phases on the shortest and middle baselines 
in the first visibility lobe are zero. This results in a predicted 
CP of k for the CO band head and leads to the jumps in the CP 
as observed. Because the model is spherically symmetric, the 
observed non-zero/non-7r CPs and non-zero DPs are not repro- 
duce^ However, the overall agreement in the spectrum, visibil- 
ities, and CP suggests that our model represents the approximate 
picture of the outer atmosphere of BK Vir. 

It is impossible to uniquely constrain the detected inhomo- 
geneous structure from our observations. However, we estimated 
the flux contribution of asymmetric features from the observed 
CPs as follows. We added a Gaussian-shaped spot to the best-fit 
model image in the CO lines and computed the visibility and CP 
at the observed baselines. We changed the size and intensity of 
the Gaussian-shaped spot so that the observed non-zero/non-7r 
CPs (40-100°) can be reproduced without affecting the fit to the 
observed visibilities significantly. This experiment suggests that 
the flux contribution of the asymmetric feature is 3-5% of the to- 
tal flux. These values agree with those estimated for other AGB 
stars. Ragland et al. (2006) derived a flux contribution of 3% 
from a geometrical modeling (uniform disk + unresolved spot) 
of the asymmetric features detected in a sample of AGB stars in 
the //-broadband. Wittkowski et al. (1201 11 1 also suggest a flux 
contribution of up to ~3% in the 2 fim H2O band for the Mira 
star R Cnc. 

The observed spectrum is nearly equally well reproduced by 
the MARCS-only model and MARCS +MOLsphere model, de- 
spite the clear disagreement between the observed visibilities 
and MARCS-only model. The reason is that the additional ab- 
sorption due to the CO layers above the MARCS model is filled 



6 The model DPs for the shortest and middle baselines are zero, and 
the model DP for the longest baseline is the same as the model CP with 
its sign flipped. 
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Table 2. Parameters of the MARCS +MOLsphere model for BK Vir 

Parameter Searched range Solution 

CO column density of the inner layer: AU (cm- 2 ) 10 21 ,2x 10 21 ,5x 10 21 , 10 22 , 10 2J (l-2)xl0 22 

CO column density of the outer layer: N outa (cm" 2 ) 5 x 10 18 , 10 19 ,2x 10", 5 x 10 19 , 5 x 10 20 10 19 -10 20 

Temperature of the inner layer: T mna (K) 1800 ... 2300 (Ar mncr = 100 K) 1900-2100 

Temperature of the outer layer: r outer (K) 1200 ... 2300 (AT 0Uter = 100 K) 1500-2100 

Radius the inner layer: Ri unt:T (R+) 1.2,1.25,1.3 1.2-1.25 

Radius the outer layer: R oum (R*) 1.5, 3.5 (ARgmm = 0-5 R*) 2.5-3.0 



0.5 



i 0.5 
> 


0.6 

5 0.4 

> 0.2 


0.15 
= 0.1 

3 

>0.05 



CU 

3200 - 



a 100 

0h 



a. 






/VW| 

3 : 


_b. Baseline = 15.9m _ 


o. B 


i i i i i i i i i i i i i i i i 

aseline = 31.8m 


_d. Baseline=47.8m 


— i — j — i — i — i — i — 




■ e. 


hil 




-\ 1 1 1 1 : 







2.295 



2.3 

A lab (^ m ) 



2.305 



Fig. 7. Comparison between the observed interferometric data 
and the MARCS +MOLsphere model shown in the same manner 
as Fig. |5] 



in by the extended emission of the same CO layers in the spa- 
tially unresolved spectrum. Therefore, good agreement between 
observation and model for spatially unresolved spectra of strong 
molecular features can be misleading, because these features can 
include the contribution from the extended outer atmosphere. 

5. Discussion 

Figure|9]shows the temperature distribution of the MARCS pho- 
tospheric model as a function of radius, together with the lo- 
cation and the temperature of the extended CO layers derived 
from our modeling. The temperature of the inner CO layer is no- 
ticeably higher than that of the uppermost layer of the MARCS 
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Fig. 8. MARCS +MOLsphere model for BK Vir shown in the 
same manner as Fig. [6] 



photospheric model, although the height is approximately the 
same. The temperature of the outer CO layer is also higher than 
or equal to that of the top of the photosphere. This suggests a 
temperature inversion in the MOLsphere of BK Vir. 

Given that our spectroscopic calibration method without a 
usual spectroscopic standard star gives the reliable CO line spec- 
tra of BK Vir as described in Sect. [2] the uncertainty in the spec- 
tra are unlikely to affect the result of our modeling. 

We checked whether this is merely a result of approximat- 
ing the temperature distribution in the outer atmosphere by two 
discrete layers. We set the temperatures of the CO layers to be 
below 1500 K, which is the temperature of the uppermost layer 
of the MARCS photospheric model, and attempted to reproduce 
the observed data. However, it turned out that such models can- 
not provide a reasonable fit to the observed spectrum and visibil- 
ities simultaneously. Therefore, the temperature inversion sug- 
gested in Fig. |9]is unlikely to be a result of the simplification we 
adopted for the structure of the outer atmosphere. 
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We also checked whether or not the assumption of local 
thermodynamical equilibrium (LTE) used in our model is valid 
for the CO layers. As described in Ohnaka (2004), we com- 
pared the collisional and radiative de-excitation rates. The col- 
lisional de-excitation rate C u \ is estimated as N cr ul u re i, where 
N is the number density of the primary collision partner, <r u i is 
the cross section, which is approximated with the geometrical 
cross section, and v te \ is the relative velocity between the col- 
lision partner and CO molecules. The CO column density of 
10 21 cur 2 derived for the inner layer and the layer's geomet- 
rical thickness of 0. 1 R* translates into a CO number density 
of 7 x 10 9 ctrT 3 . In chemical equilibrium at the temperature of 
2000 K derived for the inner layer, number densities of H atoms 
and H2 molecules of 1.2 x 10 13 and 1.4 x 10 13 cm -3 are nec- 
essary to obtain the above CO number density. We take both H 
atoms and H2 molecules as the primary collision partner of CO 
molecules, which leads to N - 2.6 x 10 13 ctrT 3 . With a geomet- 
rical cross section cr u i of 10~ 15 crrT 2 and a relative velocity u re i 
of 5 km s _1 assumed, we obtain C u i = 1.3 X 10 4 s _1 for the inner 
CO layer. In the case of the outer CO layer, the lower limit of 
the CO column density (10 19 cm -2 ) translates into a CO number 
density of 7 . 1 x 1 6 crrT 3 . For the temperatures of 1 500-2 1 00 K 
derived for the outer layer, the number densities of H atoms and 
H 2 molecules are 5.0 x 10 9 and 1.6 x 10 10 cm" 3 (1500 K) and 
3.8 X 10 10 and 4.2 x 10 7 ctrT 3 (2100K), respectively. This leads 
to N (H + H 2 ) = (2^1)xl0 10 cm" 3 . With the above cr ul and y rei , 
we obtain C u \ = 10s _1 for the outer CO layer. These collisional 
de-excitation rates are much larger than radiative de-excitation 
rates A u \ of 0.6-1 s found for the CO lines analyzed in the 
present work, which means the validity of LTE for our model- 
ing. Therefore, it is unlikely that the high temperatures derived 
for the CO layers are just a result of non-LTE effects. 

The temperature inversion implies the operation of some 
heating mechanism in the outer atmosphere within ~3 R*. While 
the outwardly propagating shocks generated by the stellar pulsa- 
tion can produce temperature jumps in Mira stars (e.g., Richter et 
al. 120031 . it is not clear whether this is the case for normal M gi- 
ants, where the variability amplitudes are much smaller. Despite 
the remarkable difference in the variability amplitude, BK Vir 



shows the temperature inversion in the molecular component of 
the outer atmosphere. This implies that the stellar pulsation may 
not be the cause of the heating. 

Some kind of heating is also required to explain the chro- 
mosphere common among M giants, although no UV or Ha ob- 
servations of BK Vir are available in the literature. However, 
the (electron) temperature of the chromosphere inferred from 
the UV emission lines in M giants is much higher (~6000 — 
10000 K, Eaton & Johnson flMgl Luttermoser et al. [T994| than 
the 1900-2100 K derived for the inner CO layer. Therefore, al- 
though the inner CO layer is warmer than the uppermost layer of 
the photosphere, the CO gas in the extended outer atmosphere 
represents a distinct component other than the chromosphere. 

Still, the process responsible for the chromospheric heating 
may lead to the formation of both of the hot plasma and warm 
molecular gas in the outer atmosphere. For example, the mag- 
netohydrodynamical (MHD) simulation of stellar winds in red 
giants by Suzuki (2007) shows the formation of transient hot 
gas bubbles embedded in cool gas outflows, although the simu- 
lation is still limited to stars with relatively high effective tem- 
peratures (^3900 K) and much less luminous stars ($1000 Lq) 
compared to BK Vir. However, as Airapetian et al. (2010) point 
out, this simulation assumes fully-ionized plasma, although the 
gas around K-M giants is only weakly ionized. The MHD sim- 
ulation of Airapetian et al. (f2010) reproduces the observed ter- 
minal wind velocities and mass-loss rates of K giants, but the 
agreement between the modeled and observed radial velocity 
fields is not yet satisfactory (comparison between the model of 
Airapetian et al. l2010l and observed data summarized in Crowley 
et al. 120091 is shown in Harper 1201 II ). Since the simulation of 
Airapetian et al. (120101 ) assumes a nearly iso-thermal chromo- 
sphere, it does not include the cool, molecular component. The 
extension of these MHD simulations to cooler and more lumi- 
nous stars including the modeling of thermal structures would 
be useful for understanding the origin of the MOLsphere in nor- 
mal red giants. 

The presence of the dense CO layers above the photosphere 
suggests that there can also be other molecular species in the 
MOLsphere. For example, TiO is one of such molecules, and 
the contribution from the MOLsphere can significantly affect the 
strength of the TiO bands observed in the optical. Hron et al. 
(2010) show that the MOLsphere models for the red supergiant 
Betelgeuse available in the literature predict the TiO bands to be 
too strong compared to the observed ones. This is because of the 
strong contribution from the TiO molecules in the MOLsphere. 
As possible reasons, they suggest the clumpy or patchy nature of 
the MOLsphere with a small filling factor and/or the crude nature 
of the current MOLsphere models. In addition, we suggest that 
the TiO bands may form by scattering in the MOLsphere. The 
line formation by scattering for electronic transitions (such as the 
TiO bands) is suggested by Hinkle & Lambert ( 1975 ). The con- 
struction of the MOLsphere model with this process taken into 
account is beyond the scope of this paper, but should be pursued 
to obtain a more comprehensive picture of the MOLsphere. 

6. Concluding remarks 

We have spatially resolved the normal M7 giant BK Vir in 
the individual CO first overtone fines with VLTI/ AMBER. The 
uniform-disk diameter shows clear increases by 12-31% in the 
CO lines compared to the continuum. We detected non-zero/non- 
7i CPs and non-zero DPs in the CO lines, which reveal an asym- 
metry in the upper photosphere and outer atmosphere where the 
CO lines form. The visibilities, CPs, and DPs observed in the in- 
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dividual CO lines are symmetric with respect to the line center. 
This is in marked contrast with the red supergiant Betelgeuse, in 
which the vigorous inhomogeneous gas motions manifest them- 
selves as the asymmetry in the visibilities in the blue and red 
wing of the lines. The absence of such asymmetry indicates that 
the amplitude of the inhomogeneous velocity field in BK Vir is 
much smaller (S=5 km s _1 ) than in Betelgeuse. 

While the MARCS photospheric model can well reproduce 
the observed CO line spectrum of BK Vir, the angular size in the 
CO lines predicted by the model is too small. This indicates the 
presence of more extended CO layers than predicted by the pho- 
tospheric model. Our model with two CO layers added above the 
MARCS model can explain the observed spectrum, visibilities, 
and CPs reasonably well. The derived temperatures of the two 
extended CO layers, 1900-2100 K and 1500-2100 K, are higher 
than or equal to that of the uppermost layer of the photosphere. 
This result suggests the operation of some heating mechanism in 
the outer atmosphere within ~3 

The present work illustrates that AMBER observations with 
high-spectral resolution are effective for constraining the phys- 
ical properties of the outer atmosphere. Given that normal red 
giants with small variability amplitudes are not yet intensively 
studied by IR spectro-interferometry, it is indispensable to carry 
out a systematic survey of a sample of these stars with ranges of 
temperatures and luminosities. This will shed new light on the 
mass-loss mechanism not only in normal red giants but also in 
Mira stars. 
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